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Abstract. 

A recently developed wavelet based approach is employed to characterize the scaling 
behavior of spectral fluctuations of random matrix ensembles, as well as complex 
atomic systems. Our study clearly reveals anti-persistent behavior and supports 
the Fourier power spectral analysis. It also finds evidence for multi- fractal nature 
in the atomic spectra. The multi-resolution and localization nature of the discrete 
wavelets ideally characterizes the fluctuations in these time series, some of which are 
not stationary. 

PACS numbers: 05.40.-a,05.45.Mt,05.45.Tp,32.30.-r 

Recently, a fundamentally different approach to the study of random matrix 
ensembles has been taken [HI21EII1I- The cumulative of the unfolded spectra for a variety 
of ensembles has been treated as a time series, on which a Fourier spectral analysis has 
been carried out. Very interestingly, this reveals a ubiquitous 1/f power law behavior in 
the Fourier domain, which has been taken as the defining characteristic of these random 
matrix ensembles. In this context, a number of tools, earlier employed for the analysis 
of various time series showing self-similar behavior, can be used for the characterization 
of the above data. The well-known methods like rescaled range analysis structure 
function method [Bj, wavelet transform modulus maxima [7j etc., although powerful, 
have been found to be wanting, when applied to non-stationary time series. Recently, 
we have developed a discrete Daubechies wavelets |B] based approach [01 ED], which 
is well suited to remove local trends in any time series and faithfully characterize the 
fluctuations. It is complimentary to the well studied multi-fractal detrended fluctuation 
analysis (MFDFA) [TTJ[T2J[Tni[IU[I3[Tn], which uses appropriate polynomial fits in local 
windows for extracting fluctuations. It is worth mentioning that, the multi- resolution 
ability of the wavelets naturally makes them ideal for the analysis of time series showing 
self-similar behavior [T7] . 

This letter studies the scaling and correlation behavior and multi-scaling properties 
of energy fluctuations in random matrix ensembles and atomic systems. The primary 
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motivation for this comes from the aforementioned Fourier spectral analysis, which is 
known to yield the Hurst exponent H, with a = 2H + 1, a being the exponent of the 
power law decay. Since, H is related to the second moment, it is quite natural to ask 
the nature of the higher moments, which determine the multi-fractal characteristic of 
the time series. For this purpose, we make use of a discrete wavelet based approach, 
developed by us earlier [UJ EH . In this approach, after removal of the local trend over a 
given window, the fluctuations are extracted. The wavelets naturally provide a number 
of windows of different sizes, to extract fluctuations at various scales. Using these, 
the fluctuation function [TH] is then calculated, which yields the mono or multi-fractal 
nature of the time series, when studied in a log-log plot. We further test the self similar 
behavior of the above time series through MFDFA for completeness. Our results on 
Hurst exponent matches with Fourier spectral analysis. 

It is found that certain atomic systems exhibit multi-fractal behaviour. It is 
observed that removing the eigenvalues of the localized states, which are quantum 
analogue of classical unstable orbits, from the energy level data does not alter the 
nature of the self-similar behavior, although it affects the Hurst exponent value. The 
identification and quantification of localized states are discussed in Ref. [TH]. Study of 
multi-fractal behaviour has been carried out earlier in the context of various atoms and 
ions. Methods like, box-counting [20] and correlation sum have been employed [2~H I22j. 
Keeping in mind the possible non-stationary nature of the data, as well as efficacy of 
the wavelet based approach for extracting fluctuations, we have used the same here for 
studying the self-similar nature of the fluctuations. 

We give a brief description of the wavelets before proceeding to the analysis of time 
series related to various random matrix ensembles [2HJ an d atomic level data [2l] . Results 
and discussions pertaining to the above analysis are then presented. One observes mono- 
fractal behavior for the time series associated with random matrix ensembles. Two 
atomic energy level time series show multi-fractal behavior and another reveals mono- 
fractal character. We summarize and conclude in the final section, after pointing out 
several directions for future work. 

The discrete wavelets provide complete and orthonormal basis functions, starting 
from the father wavelet (f)(t) (scaling function) and mother wavelet ip(t) [8.21]|8| 125]. 
These functions necessarily satisfy, J <p(t)dt = A and J ip(t)dt = 0, where A is constant. 
Scaling and translation of wavelets lead to iftj^ = 2^ 2 ip(2H — /c), which obey the 
orthogonality conditions: (0j,fc|'0j,fc) = and (V^ifclVy.fc) = bj,j'bk,k'- Any signal belonging 
to L 2 can be written in the form, 

oo oo oo 

/(*)= Yl Ck ^(*) + Yl ' (*) 

k=— oo k=— oo j=0 

where c' k s are the low-pass coefficients and d'^ k s are the high-pass coefficients. The 
Daubechies family of wavelets are made to satisfy vanishing moment conditions: 
J dt t m ipj^{t) = 0. This makes them ideal to isolate polynomial trends from fluctuations. 
We make use of the discrete wavelets from Daubechies family for our analysis of energy 
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level fluctuations. 

Consider E^, i = l,2,3...,n + 1, the discrete energy levels of random matrix 
ensembles or atomic systems, represented as a discrete time series in the form Ei as 
shown in Figs. 1 and 2. The trend and fluctuations from the integrated level density 
are separated by the higher order polynomial fitting: N(E) = N t (E) + Nf(E). The 
unfolded energy spectra E^ is obtained through the transformation N t (Ej) = E^. The 
fluctuations of the energy level spacings are, = E i+1 — E iy i = 1,2, 3,..., n, where 
< ei >= 1. The time series of energy level fluctuations is given by, 

m 

5 m = ^J(e— < ei >); m=l,...,n. (2) 

i=i 

The power spectrum of the times series 5 m is calculated for various random 
matrix ensembles and atomic energy levels of Nd, Pm and Sm. It should be noted 
that, these have 60, 61 and 62 electrons respectively. The corresponding number of 
active valence electrons are 6, 7 and 8; the complexity increases as the number of 
active valence electrons increases. The random matrix ensembles studied are Gaussian 
orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE), Gaussian diagonal 
ensemble (GDE) and Gaussian symplectic ensemble (GSE). The cumulative of the 
unfolded time series are shown in Fig.l and Fig. 2, for the atomic systems and random 
matrix ensembles respectively. Regarding the details of the atomic system, interested 
readers are referred to Ref.j21] 

Fourier spectral analysis of all the data sets yields, S(k) ~ with a ~ 1. 
This is in agreement with the results of Relano et. al.„ which found that the energy 
spectra of quantum systems exhibiting classical chaos are characterized by 1/ f behavior. 
We carried out the spectral analysis using fast Fourier transform [26j, and ensemble 
averaging was employed in finding the exact slope a. 

We now proceed to the study of the scaling properties of these time series through 
wavelet based fluctuation analysis. We make use of wavelets from Daubechies family for 
characterization, since these naturally remove polynomial trends from data sets. The 
fractal nature of the time series is revealed through the study of the fluctuation function. 
In discrete wavelet transform it is well-known that, a given signal belonging to L 2 space 
can be represented in a nested vector space spanned by the scaling functions alone. This 
basic requirement of multi- resolution analysis can be formally written as [?], 

... C z/-2 C V-i C z/_o Ci^iC u 2 ... C L 2 , (3) 

with z/_oo = and = L 2 . This provides a successive approximation of a given 
signal in terms of low-pass or approximation coefficients. It is clear that, the space that 
contains high resolution signals will also contain signals of lower resolution. The signal 
or time series can be approximated at a level of ones choice, for use in finding the local 
trend over a desired window. The fluctuations can then be obtained by subtracting 
the above trend from the signal. We have followed this approach for extracting the 
fluctuations, by elimination of local polynomial trends through the Daubechies wavelets. 
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We compute the fluctuation function in order to ascertain the self-similar nature of 
the time series. The q th order fluctuation function, F q (s) is obtained by squaring and 
averaging fluctuations over all segments: 



s b=i 

Here 'g' is the order of moments that takes any real values. The above procedure 
is repeated for variable window sizes for different value of q (except q = 0). The scaling 
behavior is obtained by analyzing the fluctuation function, 

F g (s) ~ s h ^\ (5) 

in a logarithmic scale for each value of q. If the order q = logarithmic averaging has 
to be employed to find the fluctuation function: 

F (s) = exp{— -^/n[F 2 (6, S )]}. (6) 

s 6=1 

As is well-known, if the time series is mono-fractal, the h(q) values are independent 
of q. For multifractal time series, h(q) values depend on q. The correlation behavior is 
characterized from the Hurst exponent (H = h(q = 2)), which varies from < H < 1. 
For long range correlation, H > 0.5, H = 0.5 for uncorrelated and H < 0.5 for long 
range anti-correlated time series. We refer the interested readers to j2l QDj for the details 
of this approach. The power law manifests itself as a straight line in the log-log plot of 
F g (s) versus s for each value of q: 

F q (s) ~ s h ^\ (7) 

For mono-fractal time series, h(q) is constant for all q, whereas for multi-fractal time 
series h(q) shows non-linear dependence for all q. Here q varies from —10 to +10. We 
have used Db-8 wavelet for capturing fluctuations. The well known Hurst exponent 
H equals to h(q = 2), which is related to the power spectral analysis by the relation 
a = 2H + 1. Since, the values of the fluctuations are very small, in order to study 
the same, we integrated the time series by subtracting mean. Through this double 
integrated time series, the obtained Hurst exponent is H\ = H + 1. 

The Hurst exponent calculated from Fourier analysis compares well with the wavelet 
based fluctuation analysis for correlation behavior. We have also corroborated our 
results through MFDFA. 

The self-similar behavior of a variety of random matrix ensembles and atomic 
level data have been explored through both discrete wavelets and MFDFA. The scaling 
behavior corroborates the findings of Fourier analysis. We found mono-fractal behavior 
for random matrix ensemble time series for which computed Hurst exponents agreed 
with the results of Relano et. al. Very interestingly, for one atomic level data, we 
observed multifractal behavior for Pm and Nd systems. Sm showed scaling behaviour. 
Removal of the eigenvalues of the localized states made the Sm system monofractal 
with strong persistence. These results are shown in Table. 1. In the last column Sm^ 
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Figure 1. Time series of extracted fluctuations from the unfolded energy spectra of 
[a] Pm, [b] Nd and [c]and [d] Sm atomic energy levels. 
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Figure 2. Time series of extracted fluctuations from the unfolded energy spectra of 
random matrix ensembles [a] GDE, [b] GSE, [c] GOE and [d] GUE. 



represents the atomic level data for Sm when the eigenvalues of the localized states have 
been removed. It is worth noting that, among the three atoms, Sm has the strongest 
configuration mixing as it has the largest number of active valence electrons. This 
indicates that the lack of sufficient mixing and effect of localized levels may influence 
the multifractal nature of the atomic systems. 

We intend to study these aspects carefully in future. The nature of correlations in 
other ensembles which appear in various physical problems also needs investigations. 
These include embedded Random Matrix ensembles relevant for finite interacting 
particle systems [27| . 

We are thankful to Profs. J. C. Parikh and V.K.B. Kota for many useful discussions. 
We acknowledge extensive discussions with Drs. M.S. Santhanam and D. Angom, to 
whom we are also thankful for providing us with atomic level data. 
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Table 1. The h(q) values for different values of q obtained from wavelet based 
fluctuation analysis. Here h(q = 2) = H , is the Hurst scaling exponent. 
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